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Abstract 

00 

An explicit formula is derived for the electromagnetic (EM) field scattered by one 
small impedance particle D of an arbitrary shape. If a is the characteristic size of 
the particle, A is the wavelength, a « A and £ is the boundary impedance of D, 
[A, [. E , N]\ = £[IV, H] on S, where S is the surface of the particle, N is the unit outer 
normal to S, and E, H is the EM field, then the scattered field is E sc = \S7g(x, x\),Q]- 
Here g(x,y) = 4 / 1 / / , k is the wave number, x\ £ D is an arbitrary point, and 

Q = — V/rV x Eq, where Eq is the incident field, |5| is the area of S, oj is the 
frequency, g is the magnetic permeability of the space exterior to D , and r is a tensor 
which is calculated explicitly. The scattered field is 0(|C|a 2 ) >> 0(a 3 ) as a —>• 0 
when A is fixed and ( does not depend on a. Thus, \E SC \ is much larger than the 
classical value 0(a 3 ) for the field scattered by a small particle. It is proved that the 
effective field in the medium, in which many small particles are embedded, has a limit 
as a —> 0 and the number M = M(a ) of the particles tends to 00 at a suitable rate. 
This limit solves a linear integral equation. The refraction coefficient of the limiting 
medium is calculated analytically. This yields a recipe for creating materials with a 
desired refraction coefficient. 


1 Introduction 

Let D C R 3 be a bounded domain with a connected smooth boundary S, D' := R 3 \ D, 
k 2 = const > 0 is the wave number, 00 > 0 is frequency, the boundary impedance ( = const, 
Re ( > 0, e > 0 and g > 0 are dielectric and magnetic constants, e' = e + i^, a = const > 0 
is the conductivity of D ', x G D', r = |x|, N is the unit normal to S pointing into D'. 
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Let us assume that the electrical held E = E 0 + e, where E 0 is the incident held and e 
is the scattered held. Then e solves the problem: 

V x e = iujfih , V x h = — icoe'h in D\ (1) 

r(e r — ike) — o( 1), r —> oo, (2) 

[N, [e, IV]] - V x e] = -/ on S. (3) 

Here / = [N, [E 0 , iV]] — ^[A - , V x E 0 \ is a given smooth tangential held on S, H — 
h = e' and /i are dielectric and magnetic constants of the medium IT, [A, B] = A x B 

is the cross product of two vectors, A ■ B is their scalar product. Problem (l)-(3) is the 

scattering problem for electromagnetic (EM) waves by an impedance body D of an arbitrary 
shape. This problem has been discussed in [I], where the uniqueness of its solution has 
been proved. Explicit formula for the plane EM wave scattered by a small impedance body 
(ha < 1, a is the characteristic size of this body) of an arbitrary shape is derived in [I]. 
There one can also find a solution to many-body scattering problem in the case of small 
impedance particles (bodies) of an arbitrary shape. 

A few historical remarks are in order. The theory of wave scattering by small bodies was 
originated by Rayleigh in 1871. He understood that the main term in the scattered held is 
the dipole radiation. How to calculate this radiation, in other words, how to calculate the 
induced dipole moment on a small body of an arbitrary shape, Rayleigh did not show. This 
was done nearly 100 years later by the author, see [2j. In 1908 G.Mie published a method 
for solving scattering problems for well conducting spherical particles using separation of 
variables in the spherical coordinates. His method works also for spherical impedance 
particles, but does not work for particles of an arbitrary shape for which separation of 
variables cannot be used. Smallness of the particle is not required by the Mie’s method. 

There were no analytical methods for solving EM wave scattering problems for small 
bodies of an arbitrary shape. No explicit formulas for the scattered fields were obtained by 
other authors for small bodies of an arbitrary shape. Such methods for acoustic and EM 
waves were developed in the monograph [I], 

Let us highlight the novel points in our theory: 

a) For one small impedance particle of an arbitrary shape an analytic formula for the 
scattered held is derived; this formula is asymptotically exact as a —* 0; the scattering 
amplitude is 0(C|<S'|). Assuming |*S'| = 0(a 2 ) and ( = 0(a~ K ), k G [0,1) is a constant, one 
obtains that the scattering amplitude is 0(a 2-K ) >> 0(a 3 ), where 0(a 3 ) is the value of 
the scattering amplitude in the classical theory. This is a new physical phenomenon. If the 
particle is perfectly conducting, the corresponding scattering amplitude is 0(a 3 ) for small 
a. 

b) In the case of many small impedance particles a method for solving the EM wave 
scattering problem is developed (see formulas (10)-(12)) and the limiting integral equation 
(see formula (13)) is derived for the limiting effective held in the medium where very many 
small impedance particles are embedded. 

c) These results are used to formulate a method for creating materials with a desired 
refraction coefficient, see Section 5. 

These results do not intersect with the results published by other authors. 
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2 Formula for the solution of the EM wave scattering 
by one small body 

Our first result is the following explicit formula for the scattered held: 

e ik\x-y\ 

4tt\x — y\ J ^ 

which is asymptotically, as a —* 0, exact. The quantity Q in formula (4) is given by the 
following formula: 

Q = -^rV x E 0 , (5) 

zup 

where E 0 is the incident held, the tensor r is defined as follows: 

t = I — b, b= (b jm ) = j-^r J N j (s)N m (s)ds, (6) 

and ds is an element of the surface area. In formula (5) the terms of the higher order of 
smallness, as a —* 0, are neglected. 

When we write a —> 0, it means that A is hxed: the physical meaning has the ratio a/ A. 
Formulas (4)-(6) are proved in [I]. In their derivation a new representation of the 
scattered held e is used: 

e = Vx / g(x,t)J(t)dt , / J(t)dt:—Q, 

Js Js 

where J is a tangential to S held. One can hnd e of this form if k 2 is not a Dirichlet eigen¬ 
value of the Laplacian in I). If k 2 is an eigenvalue of the Laplacian in I), then one can use 

g p (x,y) in place of g(x,y). Here g p (x,y ) is the Green’s function of the Helmholtz operator 
in the exterior of 5(0, p), the ball centered at the origin and of radius p. The g p (x,y ) 
satishes the Dirichlet condition on the boundary of 5(0, p) and the radiation condition at 
inhnity. The small number p > 0 is chosen so that k 2 is not a Dirichlet eigenvalue of the 
Laplacian in D\B(0, p). Such a ball can always be found and p > 0 can be chosen as small 
as one wishes. 

Let us discuss formulas (4)-(6). The choice of X\ G D in formula (4) is not important 
since a is small. If D is centrally symmetric body, then one can take as X\ its symmetry 
center. Formula (5) shows that the scattering amplitude A, A = Oddis'!) = 0(a 2 ) if ( 
does not depend on a. The classical dependence of the scattering amplitude on a is 0(a 3 ), 
which is much smaller than 0(a 2 ) as a —> 0. 

This fact might be of physical interest in applications. 

Formula (6) gives explicitly the dependence of the scattered held on the shape of the 
small body. 

The main physical (and mathematical) idea, used in (IJ for the derivation of the above 
results is the reduction of the scattering problem for a small particle to finding just one 
quantity Q rather than a boundary function. In the next Section the many-body scattering 
problem is discussed. 


e=[Vg(x,x 1 ),Q\, g(x,y ) 
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3 EM wave scattering by many small impedance par¬ 
ticles 

Let us formulate our basic results for EM wave scattering by many small impedance parti¬ 
cles D m , 1 < m < M , distributed in a bounded domain fb Let x rn G D m be some points. 
Assume that the number A/"(A) of small bodies (or points x m ) in any sub-domain A C 0 
is given by the formula 

Af(A) — —— [ N(x)dx(l + o(l)), a —> 0, (7) 

or K J A 

where N(x ) > 0 is a continuous in O function, k G [0,1) is a parameter, and the boundary 
impedance of the m— th small particle is Q m := h(x m )/a K , where h(x) is a continuous 
function in O, such that R eh > 0, and x rn G D m is an arbitrary point. 

The choice ( = h(x)a~ K is not dictated by physical laws. As we mentioned earlier, the 
choice of ( is made by the experimentalist as he wishes. 

The only physical restriction on the boundary impedance is the relation Re( > 0, that 
is, Reh(x ) > 0. 

The restriction for the parameter n G [0,1) is of technical nature and is not related to 
a physical law. 

The functions N(x), h(x) and the parameter k can be chosen by the experimentalist as 
he/she wants. 

Our main physical assumption is: 


a « d « A, (8) 

where d is the minimal distance between neighboring small particles. 

Our second result is the following formula for the solution of many-body EM wave 
scattering problem: 


M 

E{x) = E 0 (x) + y^ j [Vg(x,x m ),Q m ], a0, (9) 

m= 1 

where Q m are defined by the formula: 

Qm = -^^rV x E em . (10) 

lU)pi 

Here we assumed for simplicity that the particles have the same shape, so tensor r does 
not depend on m, and E em is the effective held acting on the m —th particle. This held is 
defined by the formula: 


M 

E e m Eo(x rn fj T ^ ^ [ ^9 (*Hn, 3Cj ), Qj\• 


( 11 ) 
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Formulas (10)-(11) lead to the linear algebraic system for finding the unknown quantities 
E m • F/ em . 




M 

-—~~ [\/g(x m , Xj), rV x Ej}hja 2 ~ K , 1 < m < M, 

" 3=hj 


( 12 ) 


where Cs > 0 is the constant in the formula \S\ = c^a 2 . 

Our third result is the following integral equation for the limiting effective field in 0: 


E(x) = E 0 (x) - —V x [ g(x,y),rV x E(y)h(y)N(y)dy, (13) 

J n 

where N(y) is defined in formula (7) and the limit is taken as a —» 0. The existence of this 
limit is proved in [Tj. 

Equation (13) is equivalent to the following local differential equation: 


V x V x E(x) 


k 2 E(x) - —V x (h{x)N(x)rX7 x E(x)\ 

lUJ l~L \ / 


(14) 


as one can check by applying the operator V x Vx to equation (13), see the details in [I], 
Let us assume that r is proportional to a diagonal matrix I. This happens, for example, 
if the particles are balls of radius a, in which case r = 1 1, as one can easily verify. 

Then equation (14) takes the form: 


V x V x E(x) 


k 2 E(x) 2c s [v(/i(^)),VxL(i)] 

1 + £%; h ( x ) N ( x ) 3iuj V 1 + i^ h ( x ) N ( x ) 


(15) 


4 Physical interpretation of formula (15) 

To interpret physically formula (15), consider the Maxwell’s equations: 

V x E(x) = iuogH, V x H{x) = —ioue'H, (16) 


where /i = g(x), apply the operator Vx to the first equation and then use the second one. 
This yields: 

V x V x E(x) = K 2 E{x) + V x E(x) ], (17) 

g{x) 

where K 2 = u 2 g(x)e' = k 2 ^^ , k 2 = Lu 2 e 0 g 0 , and /U 0 ,e 0 are the parameters of the free 

space. Let n(x) := . 

^ v ' m on) 

Comparing formulas (15) and (17) one concludes that the limiting medium, obtained 
by the embedding of many small impedance particles, has the new refraction coefficient: 


n(x) 


n 0 (x) 





(18) 
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and the new magnetic permeability. 


H{x) — (l +-^-h(x)N(x)] . (19) 

V diup, / 

On formula (18) our recipe for creating materials with a desired refraction coefficient is 
based. This recipe is discussed in the next Section. 


5 A recipe for creating materials with a desired re¬ 
fraction coefficient 


Let us rewrite formula (18) as 


n(x) 


n 0 (x 


1 — ic\h{x)N [x] 


1/2 • 


( 20 ) 


where C\ := > 0 and h = h\ + ih 2 , hi > 0. The functions h(x) and N(x) > 0 are at 

the disposal of the experimentalist, as was mentioned earlier. By choosing these functions 
properly, one can get any desired refraction coefficient which has the property Im n(x) > 0. 
One rewrites formula (20) as follows: 


n(x) — n 0 (x) ^1 — icih(x)N(x)j = n 0 (x) ^1 + Cih 2 {x)N{x) — icihi(x)N(x )j 2 . (21) 


Let z be a complex number. Define z 1 ^' 2 = \z \ 1 , /2 e *^/ 2 , where 0 is the argument of z, 
0 < 0 < 27t. By choosing h 2 < 0 so that 1 + Cih 2 (x)N(x) is small and choosing hi > 0 
suitably, one can make |n(a;)| to be any desired non-negative function. The argument 0 of 
the expression 

1 + cih- 2 {x)N {x) — ic\hi{x)N(x) 

can be made arbitrary by choosing hi > 0 and h 2 G (— oo, oo ) suitably. 

For example, assume that n 0 > 0, hi > 0 is small and 1 + Cih 2 {x)N(x) > 0. Then the 
0 = 27r — 25, where 5 is arbitrarily small, 0/2 = 7 r — 5 and n(x) = |n(x)|e _ * < ' 7r_<5 / Thus, Re 
n(x) < 0, Im n(x) > 0, and Im n(x) = \n(x)\ sin 5 can be made as small as one wishes if 5 
is sufficiently small. Therefore, one gets a material with negative refraction and negligible 
losses. 

Similarly, using formula (19), one can change magnetic permeability in a desired direc¬ 
tion by embedding in a given medium many small impedance particles with the suitable 
boundary impedances. 
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